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Abstract 

Nucleon electromagnetic form factors are studied in the cloudy bag model 
(CBM) with center-of-mass and recoil corrections. This is the first presenta¬ 
tion of a full set of nucleon form factors using the CBM. The center of mass 
motion is eliminated via several different momentum projection techniques 
and the results are compared. It is found that the shapes of these form fac¬ 
tors are significantly improved with respect to the experimental data if the 
Lorentz contraction of the internal structure of the baryon is also appropri¬ 
ately taken into account. 


I. INTRODUCTION 

Form factors characterise the internal structure of subatomic particles and, in partic¬ 
ular, electromagnetic probes of hadrons provide important information on the underlying 
quark and gluon degrees of freedom. In the nonperturbative regime (i.e., at low momentum 
transfer), QCD-motivated, effective hadronic models continue to play an important role in 
analyzing and understanding a wealth of experimental data. The MIT bag model |l[ was 
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an early attempt to include the key features of confinement and asymptotic freedom in a 
quark based model of hadronic structure. The cloudy bag model (CBM) [S] improves on the 
MIT bag model significantly by introducing an elementary pion held coupled to the quarks 
inside the bag such that chiral symmetry is restored. The introduction of the pion held not 
only improves the static nucleon properties, but also provides a convenient connection to 
the study of conventional intermediate energy physics such as nN and NN scattering. 

There are many calculations of the nucleon electromagnetic form factors within different 
hadronic models. Indeed the understanding of these form factors is extremely important in 
any effective theory or model of the strong interaction. However, there is, to our knowledge, 
no truly satisfactory means of forming fully Lorentz covariant momentum eigenstates from 
any static model. In this work we suggest an improved treatment (a hybrid method of 
Galilean momentum projection combined with an appropriate Lorentz contraction) for a 
model which has been widely used for many years — the CBM, and bring it to larger 
momentum transfers. The present results not only remind us of the effectiveness of chiral 
quark models at moderate momentum transfer, but also remind us that they can serve as 
an essential first step in the investigation of the electromagnetic interaction in quark based 
nuclear models, in particular, the electromagnetic interaction in the quark-meson coupling 
(QMC) model ||. 

In the CBM, as in the MIT bag model, quarks are independent particles confined in a 
rigid spherical well. The bag model wave function for a baryon is a direct product of indi¬ 
vidual quark wave functions, analogous to nuclear shell model wave functions. A static bag 
cannot carry a definite momentum and so bag-model baryon states are not total momentum 
eigenstates, in spite of the fact that the Hamiltonian commutes with the total momentum 
operator. Matrix elements evaluated between such states contain spurious center of mass 
motion. This defect compromises some of the predictive power of the model, such that only 
observables involved in very low-momentum transfer processes (q 2 /Am 2 N <C 1) are typically 
assumed to be reliable. 

Over the years, a number of prescriptions for the correction of the center of mass motion 
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have been developed (for an overview, see for example Refs. [[|]). The diversity of approaches 
may be viewed as an indication of the uncertainty associated with this correction. In contrast 
to the nonrelativistic case, the internal motion of a composite object cannot be explicitly 
separated from the collective motion in a covariant description. For the calculation of the 
form factors, a satisfactory treatment may result from a combination of relativistic boost, 
momentum projection, and a variational procedure. Betz and Goldflam [j|] argued that a 
static soliton bag can be boosted consistently to a soliton bag moving with a finite velocity. 
However this approach is impractical for boosting the MIT bag because of the sharp surface 
which prevents the construction of a simple boost operator || . A number of nonrelativistic 
methods for the center of mass correction exist in the literature [7]M|9|. Analytic forms of 


the recoil corrections can be obtained in a relativitic harmonic oscillator quark model [10 


Unfortunately different groups do not always agree with each other and sometimes even 
result in a correction with the opposite sign. 

In this work we compare several intuitively simple momentum-projection procedures for 
the calculation of the nucleon electromagnetic form factors. The basic idea is to extract the 
momentum eigenstates from the static solutions by appropriate linear superpositions. The 
simplest prescription for this approach was proposed by Peierls and Yoccoz (PY) [[TTJ . We 
will assume that a baryon is composed of three constituents. Hence the wave function for a 
moving baryon with total momentum p is constructed as 


v bp Y (x 1 , x 2 , x 3 ; p) = Y py (p) /d 3 xe ,p ' x T(x 1 ,x 2 ,x 3 ;x), 


( 1 ) 


where Ypy(p) is a momentum dependent normalization constant. The localized state is 
simply given by a product of the three individual quark wavefunctions, 


T(xi, x 2 , x 3 ; x) = g(xi - x)g(x 2 - x)g(x 3 - x), 


( 2 ) 


where x refers to the location of the centre of the static bag, and xi, x 2 , and x 3 specify 
the positions of the three constituent quarks. With the PY wave function, the predictions 
of the static baryon properties are generally improved [Q. It reduces the r.m.s. radius, 
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increases 5 % and on the whole produces a better mass spectrum. However, it is unreliable 
for calculations of dynamic observables which involve large momentum transfers, sine the 
PY wave function does not transform appropriately under Lorentz boosts. 

A closely related method for eliminating the center of mass motion is called the Peierls- 
Thouless (PT) projection [|T^,|T3|]. There the wave function is constructed through one further 
linear superposition in terms of the PY wavefunction, 


Tp T (xi,X 2 ,x 3 ;p) = N(p) [d 3 p'w( p')e* (p p 'f Xc m T PY (xi, x 2 , x 3 ; p'). 


( 3 ) 


where x c . m . = (x x + x 2 + x 3 )/3 is the center of mass of the baryon (we assume equal mass 
quarks here). Ideally the weight function, w;(p'), should be chosen to minimize the total 


energy, but this is quite complicated to implement in practice. As in Ref. ||T3[, we make 
the choice u^p') = 1 for simplicity and convenience. Then integrations over x and p' can 
be carried out explicitly. This leads to a comparatively simple PT wave function for the 
baryon, 


T pt (x 1 ,x 2 ,x 3 ; p) = Ap T e* p ' Xc m g(xi - x c . m .)g(x 2 - x c . m .)g(x 3 - x c . r 


( 4 ) 


where Np t is determined by the requirement that it satisfy the normalization condition 

j d 3 Xi<i 3 x 2 (i 3 a; 3 ^ PT (x 1 ,x 2 ,x 3 ;p')^ P T(xi,x 2 ,x 3 ;p) = (27r) 3 5 (3) (p / - p)r/ p . (5) 

with r] p = 1 for the nonrelativistic normalization and r] p = E{p)/mN if we wish to adopt a 
standard relativistic normalization for the baryon wavefunction. 

Notice that the above methods of momentum projection act only on the center of mass 
coordinate and the individual quark wavefunctions are not affected. However, since baryons 
are composite objects, once they have nonzero momentum, their internal structure should 
be subsequently modified. For example, the bag surface is no longer spherical in the Breit 
frame, rather it should be contracted along the direction of motion. We take care of this 
effect in terms of the preseprition by Licht and Pagnamenta [|14[ . 

It should be noted that the present work is the first presentation of calculations of 
the nucleon electromagnetic form factors using the CBM, besides the obvious improvement 
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of the treatment. The outline of the paper is as follows. Firstly, we briefly review the 
electromagnetic interactions of the CBM in Sec. II. The calculation of electromagnetic form 
factors for the bare bag with momentum projection is then presented in Sec. III. In Sec. IV, 
we discuss the necessary scaling of the form factors due to the effects of Lorentz contraction. 
Pionic corrections are then given in Sec. V. The numerical results are presented and discussed 
in Sec. VI before the concluding remarks in Sec. VII. Some technical details and explicit 
proof of gauge invariance of the calculations are provided in appendix. 


II. ELECTROMAGNETIC CURRENTS IN THE CBM 


The linearized CBM Lagrangian with the pseudoscalar pion-quark coupling (up to order 
l/u ) is given by §J 


C = 


{iqYd^q ~ B)9 V - ^qq8 s 
+^(^ 7r ) 2 - l 771 ^ 2 ~ Jf^5 T ' 7 rq 5 s . 


( 6 ) 


where B is a bag constant, is the 7r decay constant, 9y is a step function (unity inside 
the bag volume and vanishing outside) and 8s is a surface delta function. In a lowest order 
perturbative treatment of the pion field, the quark wave function is not affected by the pion 
held and is simply given by the MIT bag solution JTj] 


g(r) 


\ 


g(r) 

\i(r-rf(r) } 


cj)9(R-r), 


(7) 


where </> contains the spin-isospin information for the wavefunction of the quark, a is the 
usual Pauli spin operator, and R is the spherical bag radius. For the ground state of a 
massless quark g{r ) = N s j 0 (co s r/ R), f(r) = N s ji(u s r /R), where co s = 2.0428 and iV s 2 = 
uj s /8ttR z jl( uj s ) ( u; s - 1). 

From the CBM Lagrangian given in Eq. ([|), the conserved local electromagnetic current 
can be derived using the principle of minimal coupling d^ —» d^ + iqA^, where q is the charge 
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carried by the field upon which the derivative operator acts. The total electromagnetic 
current is then 


j»( x )=j^ Q \x)+j^\x), ( 8 ) 

f {Q \x) = J2Qfeq f ( x )Yqf( x ), ( 9 ) 

/ 

j^\x) = —ie[ n\x)d^n(x) — n(x)d^n\x)}, (10) 


where qj(x ) is the quark field operator for the flavor f,Qf is its charge in units of e, and 
e = \e\ is the magnitude of the electron charge. The charged pion field operator is defined 
as 


tt(x) 


1 

71 


[7Ti (a;) + f7r 2 (x)], 


( 11 ) 


which either destroys a negatively charged pion or creates a positively charged one. 

The physical baryon state is then a dressed bag, consisting of a superposition of a bare 
bag and a bag with a pion cloud. Algebraically, it has the form 


\A) = y[zf[ 1 + (m A -Hq- AHjA^Hj] |A 0 ), 


where Z A is the bare baryon probability in the physical baryon states, 




i + E 


7 , 


AB \ 2 

0 \ 


dk k 4 u 2 (kR ) 


m-n 


12 tt 2 Jo uJ k (m A -m B - u k ) 2 


( 12 ) 


( 13 ) 


where A is a projection operator which annihilates all the components of A) without at 
least one pion, and Hj is the interaction Hamiltonian which describes the process of emission 
and absorption of pions. We follow the traditional CBM treatments and consider only states 
with at most one pion. The matrix elements of Hj between the bare baryon states and their 
properties are then given by JT5[ 


v$f(k) = (AolHjln^Bo) = 


wi B {k) = (Aon j (k)\H I \B 0 ) = v BA {k) = -v AB {k) = v AB (-k), 


7o u(kR ) S^riSBmsAJz* t\r<tBnt A (}* Tt \ 

k ) C T B lT A {tn tj), (14) 
( 15 ) 


BA/ 


,AB/ 


,AB / 


where the pion has momentum k and isospin projection j. Note also that f$ B is the reduced 
matrix element for the tcBq —> Ao transition vertex, u(kR ) = 3ji(kR)/kR, 04 - = \Jk 2 + m 2 , 
and s m and t n are spherical unit vectors for spin and isospin, respectively. 
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III. MOMENTUM PROJECTION CALCULATIONS FOR A BARE BAG 


It is customary to define the nucleon electric (Ge) and magnetic (Gm) form factors in 
the Breit frame by 

W'(§U°(0) |jv,(-|)> = xl' X, Ge(9 2 ), (16) 

<W>(f)|T(0)|JV.(-|)> = xt (17) 

where Xs and xl' are Pauli spinors for the initial and final nucleons, q is the Breit-frame 
three momentum transfer, i.e., q 2 = q$ — q 2 = —q 2 = — Q 2 . We choose q to define the z-axis. 
The major advantage of the Breit frame is that Ge and Gm are explicitly decoupled, and 
can be determined respectively by the time and space components of the electromagnetic 
current operator JU 

In the definition above [i.e., Eqs. ( |T6| ) and ([L7])], both initial and final states are physical 
states. Using Eqs. (||) and (|T2D, the total electromagnetic form factors can be expressed in 
terms of the three processes shown in Fig. 1. In this section we calculate the contribution 
from the bare bag only, and leave the pion loop effects to be included in a later section. 
For the three-momentum eigenstates, Eqs. ([]]) and (f|), we can proceed to calculate the 
electromagnetic form factors in a relatively straightforward way. For the PY projection, we 
obtain, 


G P e Y (Q 2 ) = I E (q 2 )/D PY (q 2 ), 

G PY (q 2 ) = I M (q 2 )/D PY (q 2 ), 


where 


W) 

W) 


D PY (q 2 ) 


poo /• 

dzz 2 Nq(z) Jd 3 rj 0 (qr) 


9+9- + 


f+f- 

r +r _ 



2m N rdzz'N+z) fd+ilM 

Jo J qr 

2 (9-f+ , 9+f-\ , r-z(g-f+ 

r {— + —) + — {— 

r oo 

J o dzz 2 N^(z)j 0 (qz/2)e, 



e, 


(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 
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0 . 


(23) 


Nq(z) = j d 3 r 


9+9- + 


f+f- 

r + r_ 




Here Dpy(q 2 ) is the momentum dependent normalization factor and Nq(z) is the overlap 
integral associated with each quark spectator. The following shorthand notation has been 
used in the above equations, 


r± = 



J 


9 ± = g(r±), f± = f(r±), 


e=e(R-r+)e(R-r.). 


(24) 

(25) 

(26) 


Similarly, for the PT projection, we obtain 

G F E T (q 2 ) = J d 3 rj 0 (qr)p(r)K(r)/D PT , ( 27 ) 

G P J{q 2 ) = J 'd 3 rj 1 (qr)g(r)f(r)K(r)/D PT , ( 28 ) 

where 

D pt = Jd 3 rp(r)K(r), ( 29 ) 

with p{r) = g 2 (r ) + / 2 (r) and K(r ) = / d 3 xp(x)p(—x — r) is the recoil function to account 

for the correlation of the two spectator quarks. 

As expected, without the momentum projection, Eqs. (|l8| , |T9l ) and Eqs. (|27|, would 
reduce to the familiar results for the static, spherical MIT bag, i.e., 

(static) (g 2 ) = I d 3 rj 0 (qr) [ g 2 {r ) + / 2 (r)], (30) 

G^c ){q 2 ) 2mN r d 3 r JiM [2 g(r)f(r)}. (31) 

J q 

Note that nonrelativistic normalization of the nucleon wave functions has been used here 
[see Eq. (^)], as is appropriate for the simple (Galilean) three-momentum projections being 
used here. 
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IV. CORRECTIONS FROM THE LORENTZ CONTRACTION 


A complete solution of a covariant many-body problem is extremely difficult. There is a 
substantial body of literature which uses light-cone dynamics [|RJ for the constituent quarks. 
With a few parameters this approach can reproduce experimental data over quite a large 
momentum transfer range. For the bag model, there is no Lorentz covariant solution for an 
extended quantum object in more than two dimensions Q- Thus we use a semiclassical 
prescription here. 

As mentioned in the introduction, the spherical bag is expected to undergo a Lorentz 
contraction along the direction of motion once it acquires a momentum. An intuitive pre¬ 
scription by Licht and Pagnamenta [bf] suggested that, in the preferred Breit frame, the 
interaction of the individual constituents of a cluster with the projectile may be regarded 
as instantaneous to a good approximation. Relativistic form factors can be simply derived 
from the corresponding nonrelativistic ones by a simple substitution rule. In the case of 
the bag model, once the spurious center of mass motion is subtracted using the PY or PT 
procedure, it is natural to rescale the quark internal coordinates as well, i.e., 


'S(x 1 ,x 2 ,x 3 ]0) <1U J , : pj 


(32) 


where the quark coordinates x\ for the moving bag are related to the Xi by a Lorentz 
transformation. Without loss of generality, we again choose the photon momentum q along 
the ^-direction. Then in the Breit frame, the quark displacements must contract in the 
^-direction while they remain unchanged in the x and y directions. Thus we have 


Za = 


m N 


E 

m N 


Zi , 


<Px\<Px' 2 = 


(33) 

(34) 


where we have assumed t = 0 for all constituents, i.e., the instantaneous approximation. 
Note that nix is the nucleon mass and E is the on-shell nucleon energy in the Breit frame. 
The (rriN/E) 2 factor is due to the Lorentz contraction of the coordinates of the two spectator 
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quarks along the direction of motion. As an example, the proton charge form factor in the 
PT scheme is thus 


G e(<? 2 ) = Jjld^x'te m ' 3 [ql(x' 3 - x' c . m .)S(x’ 3 )qf(x '3 - l' c .m.)] Pp(x'l - X'c.m)Pp{x r 2 - 4rr, l 
= (^jr) 2 ji 3 Xid z X 2 e-' qzim,,IE) p(x x - x c . m .)p(x 2 - x c . m .)p(-x c . m .) 

(35) 


‘^ffGttfrnllE 2 ), 


where qp is the quark wave function in the Breit frame (in a deformed bag), pp is the proba¬ 
bility density of the quark, and is the charge form factor calculated with the spherical 
static bag wave function [such as Eqs.([l8|) and (|27|) 1. In the second step of the derivation 
we have used the fact that a probability amplitude is a constant in different Lorentz frames, 
hence, the identity qp(x') = q(x) has been used as in Ref. |L4|. For the magnetic form factor, 
a similar expression can be derived, Gm(? 2 ) = (itin / E) 2 G s ^ h (q 2 m 2 N / E 2 ). Note that we used 
the fact that all three quarks have the same spatial wave function in obtaining Eq. (fU|) . 
The scaling factor in the argument is due to the coordinate change of the struck quark and 
the factor in the front, (mx/E) 2 , comes from the reduction of the integral measure of two 
spectator quarks in the Breit frame. Note that this prescription is similar, but not identi¬ 


cal, to the Lorentz contraction arguments used in the Skyrme model |18|. The difference is 
the (?71 tv/ E) 2 factor in front of G s ^ h which is absent in Ref. [|l8| . It might be argued that, 
since we use only a nonrelativistic momentum projection [i.e., r/p in Eq. @ cannot be fixed 
unambiguously], this factor is not well determined. 


V. CORRECTIONS FROM THE PERTURBATIVE PION CLOUD 


In the CBM with a bag radius above 0.7 fm the pion field is relatively weak and the pionic 
effects can be included perturbatively |19|]. As usual we assume that there is no more than 
one pion in the air. There are two processes contributing to the nucleon electromagnetic 
form factors due to the pion cloud. One involves the direct 77T7 t coupling shown in Fig. 1(c), 
and the other is the 'yqq coupling inside a pion loop, as in Fig. 1(b). 
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Fig. 1(c) actually contains three time-ordered subdiagrams, and has been evaluated in 


the CBM by Theberge and Thomas |H|. They gave 


G 


(tt) 




'((tt) 


(36) 


where the two terms correspond to two cases with different intermediate baryons (N and 
A). For completeness, we quote the individual contributions explicitly, 


g£V;JV) = 


7 


NN\ 2 


G?V; A) = 


367T 3 y m, 

1 (f 


d 3 k 


u(kR) u(k'R ) k • k' 


(N\t 3 \N), 


N A\ 2 


72tt 3 V 




UJkUJk' (u>k + Wfc/ 

tt(fcfl) u(CP) fc • $ 

(a; a at + Uk)(uAN + Uk')(uk + u>k> 


(N\t 3 \N), 


G£V;A 0 = 


g£ ) ( 5 2 ;A) = 


2mjy / / 

72tt 3 


AW\ 2 


ld 3 k 


u(kR) u(k'R) (q x A :) 5 


m„ 


(iVkUJk 1 


(N\t 3 \N), 


2m N l / 


WAX 2 


288tt 3 V 


ld 3 k 


m„ 


(cuAiv + 04 + cnfc') u(kR) u(k'R ) (g x A ;) 2 
(cUAiV + ^k){f^AN + l^k'ji^k^k'jy^k + a 4-' 


(37) 

(38) 

(39) 


(1Y|t 3 |1V), (40) 


where k' = k + q, tv bn — — m^, f NB is the renormalized 7 rAhB coupling constant, and 

r 3 is the third nucleon isospin Pauli matrix. 

Corresponding to Fig. 1(b), the transition matrix element can be written as 


(W/2) I J MQ) (0) | N(-q/2)) =Y, d 3 k 

BC,j J 

x (iV(f)| Ig^XOCp), TTj^)!ZTf |iV(-D> 

(uj bn + Wfc) i^CN + UJk) 

where p' = (g/2) + k and p = —(g/2) + k are the momenta for the intermediate baryons B 
and C. With the dynamical baryons and pion here, we have to evaluate the electromagnetic 
matrix elements for the intermediate processes in an arbitrary frame. Thus the matrix 
elements of J° might contain both G E (q 2 ) and Gm(q 2 ), as do those of J. It is convenient to 
use the identity 


(p'\f [Q) mp) = u(p') 


rFi(q 2 ) + 


ia ,llJ q v 

2mN 


F2(q 2 ) 


u(p), 


(42) 


where F x (q 2 ) = (G E (q 2 ) + r]G M {q 2 ))/( 1 + v) and F 2 {q 2 ) = (G M {q 2 ) ~ G E {q 2 ))/{ 1 + v) 
with 77 = —q 2 /Am 2 N . Both F\(g 2 ) and F 2 (q 2 ) are Lorentz scalar functions and hence can 
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be evaluated in any frame. However, it can be shown that after integrating over the loop 
momentum, k, the time (Ge) and space (Gm) components of Eq. ([41]) decouple again as 
long as the overall matrix element is evaluated in the Breit frame. The detailed expressions 
are messy and are therefore given in the appendix. 


VI. RESULTS AND DISCUSSION 


In this work, we have adopted the usual philosophy for the renormalization in the CBM, 
using the approximate relation, f AB ~ f NN ■ There are uncertain corrections on 

the bare coupling constant f^ N , such as the nonzero quark mass and the correction for 
spurious center of mass motion. Therefore, we use the renormalized coupling constant 
in our calculation, f NN ~ 3.03, which corresponds to the usual nNN coupling constant, 
fnNN - 0-081. 

It should point out that there is no unambigious way to implement strict gauge invariance 
(the Ward-Takahaski identities) for a composite particle [[hj. In this work, we ensure a 


somewhat weak requirement - electromagnetic current conservation [[H|]. The explicit proof 
is given in the appendix. Recall that Eqs. (|37|) , (|38|), (|39|), and (|40|) are results evaluated 
under a heavy baryon approximation. Ideally, Fig. 1(c) should be evaluated on the same 
footing as Fig. 1(b). Numerical calculations show that the recoil effects for intermediate 
baryon in the pion loop are negligible, and we may therefore ignore this recoil and use the 
standard static CBM results for the pionic correction. Consequently, the charge form factors 
at zero momentum transfer automatically satisfy the requirement of charge conservation, i.e., 
Ge(0 ) = Ge > \0 ) + G^^O) = ejv, where e jv is 1 for the proton and 0 for the neutron. 

The magnetic moments are simply the values of the magnetic form factors at zero mo¬ 
mentum transfer, /j = G^\ 0) + G$( 0). Note that the expression for the contribution from 
Fig. 1(b) is somewhat different from Ref. [[15]. Here, there is no Z 2 factor for Fig. 1(b) in 
consistence with the charge conservation. As a result of this choice the numerical contribu¬ 
tion from Fig. 1(b) increases by roughly 30% (— 0.73 for R = 1 frn ), bringing the total 
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nucleon magnetic moments a few percent closer to the experimental data. 

Table I gives the nucleon magnetic moments in this calculation. The center of mass 
correction reduces the static values by 5 — 10%. This is in contradiction with Ref. 0 but is 
consistent with Refs. fi|.|22|| . The bag radius dependence is significantly reduced by the pion 
cloud in the CBM, with little variation over the range R = 0.8 — 1.0 fm. The deficiency of 
the nucleon magnetic moments may be attributed to the higher order pionic corrections and 
explicit vector meson contributions. 

The difference between the two choices of normalization of the wave function [i.e., the 
factor rj p = 1 or E/m n in Eq. (||) ] is not significant with respect to the shape of the 
form factors. However it will smoothly scale these form factors. For the proton charge 
form factor, for example, the relativistic normalization raises the form factor roughly 5% at 
Q 2 = 0.5 GeV 2 and 10% at Q 2 = 1.0 GeV 2 . For clarity, we have always used 7j p = 1 in the 
following figures as previously stated. 

The characteristic effect of the center of mass correction on the charge form factor of the 
bare proton bag is illustrated in Fig. [2| with the bag radius R = 1 fm. The “dipole” refers 
to the standard dipole fit, F(Q 2 ) = 1/(1 + Q 2 / 0.71 GeV 2 ) 2 . The bare charge form factors 
calculated with the static bag usually drop too quickly. With the PY projection procedure, 
the form factor at moderate momentum transfer ( Q 2 ~ 0.5 GeV 2 ) increases nearly 100%. 
However the shape of the form factor does not change very much. It is generally too stiff 
and drops too fast which is mainly due to the sharp surface of the cavity approximation 
and lack of translational invariance of the wave function. Using the translational invariant 
PT projection procedure leads to improved behavior of the form factors. In particular, after 
including the correction arising from Lorentz contraction, the shape of the form factors is 
significantly improved. It is reassuring to see that the combination of Lorentz contraction 
and Galilean (nonrelativistic) momentum projection is less scheme dependent than the mo¬ 
mentum projection alone, —e.g., compare the pairs of curves PY and PT with PYL and 
PTL in Fig. 2. 

Fig-0 and Fig. 0 show the individual contributions to the charge form factors from the 
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quarks and the pion cloud, for a typical bag radius of R = 1 fm. For the proton charge 
form factor, it is clear that the 'yqq coupling terms dominate the form factor where the 
bare photon-bag coupling contributes nearly 75%. The correction from the 77T7T coupling 
decreases very quickly as the momentum transfer increases. A smaller bag radius will lead 
to a larger pionic contribution. For the neutron charge form factors, the contribution from 
the photon-bare bag coupling [Fig. 1(a)] vanishes due to the SU(6) structure. The severe 
cancellation between Fig. 1(b) and Fig. 1(c) results in a small but nonvanishing neutron 
charge form factor, with a negative mean square radius as one would expected simply from 
the Heisenberg Uncertainty Principle. With the bag radius R = 1 fm, we obtain the neutron 
charge r.m.s. radius of ( r 2 ) En = —0.14 fm 2 , to be compared with the experimental value of 


-0.12 fm 2 p4 . 

Figs. [5|-|8] show the bag radius dependence of the nucleon electromagnetic form factors. 


A large bag radius always leads to a softer form factor. We have used PT wave functions 
with Lorentz contraction in these calculations. The predictions show quite a reasonable 
agreement with the experimental data to much larger values of the momentum transfer 
than one has tended to expect. 


VII. SUMMARY 

We have calculated the electromagnetic form factors of the nucleon within the CBM, 
including relativistic corrections in the form of momentum projection and the Lorentz con¬ 
traction of the internal structure. Electromagnetic current conservation is ensured in this 
calculation which is performed in the Breit frame. This is the first time that a presentation 
of all the nucleon electromagnetic form factors has been made for the CBM. The two differ¬ 
ent procedures of momentum projection for the spurious center-of-mass motion give results 
which are relatively close to each other when Lorentz contraction effects are included. The 
Galilean invariant PT projection is generally a little better than the PY method in that 
it leads to a shape more closely resembling the dipole form. Including the corrections for 
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center of mass motion and Lorentz contraction, the numerical predictions are in rather good 
agreement with data in the region Q 2 < 1 GeV 2 . This is quite a remarkable result when 
one realizes the simplicity of the model. In particular, there are no explicit vector meson 
contributions and one possible future development would be to include txtt interactions. 

The authors are very grateful to I.R. Afnan and K. Tsushima for valuable discussions. 
We also thank S. Platchkov for sending the data for the neutron charge form factor. This 
work was supported by the Australian Research Council. 


APPENDIX A 


In Fig. 1(b), the intermediate 7 BC vertex can no longer be in the Breit frame. A 
strightforward evaluation gives the matrix elements 

k 2 - q 2 /4 


u {p'n u {p) — A T'N 


1 + 


(. E' + m N )(E + m N ) 
2 /O I t x „2 


+ 0(k), 


N'N ( q 2 /2 + k ■ q -q 2 /2 + k-q ■ 

u{p )-z - u(p) = -7 - —FT, -—- + 0{k), 


2m n 


2 rn,N \ E' + 


u(p')^ a u(p) = -{a x q) a N'N 


E + mN 
1 


+ 


qu u(p) = x 0)aN'N 


E' + mN E + mN 
q 2 /4 + k 2 


+ 0(k), 


1 - 


+ o{k), 


(Al) 

(A 2 ) 

(A3) 

(A4) 


2mN X1 ' 2 mN X \_ {E’ + 7tin)(E + m^v) 

where p’ = (q/ 2) + k, p = —(q/2) + k, E = (p 2 + m 2 N ) l ^ 2 1 E" = ( p r + m^) 1 / 2 , and the 
normalization constants are N = [(E + ?tin)/ 2mN] 1 ^ 2 , N' = [(E 1 + ?un)/ 2mN] 1 ^ 2 ■ The index 
a in Eqs. (|A3|) and (^) denotes a space component, the 0(k) terms in all equations refer to 
other pieces which are odd in k and will vanish after integration over the loop momentum k. 
Substituting the above matrix elements into Eq. ([H|) and performing some spin and isospin 
algebra, we obtain the nucleon electric and magnetic form factors originating from the 'yqq 
coupling [i.e., the combination of Fig. 1(a) and 1(b) with a proper normalization], 


(AS) 


G { S\q 2 ) = Z,_Gt\q 2 ) 

M 

+ Enn(q~) 

M 

+ A'aa (q 2 ) 

' 4/3' 




A 3 , 


,-i/sJ 
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Gg>(9 2 ) = Z 2 gS(9 2 ) 


‘ 1 i 


V 


-2/3 


Mnn(q 2 


V 


1/27 

-4/27 


M AA (q 2 


20/27 

-5/27 


\ 


+ M A n{q 


2x f 16v/2/27 ^ 


V 


-16\/2/27 


(A6) 


where the upper and lower coefficients refer to the proton and neutron respectively. Here 
G ( B and G® are the bare form factors calculated in Sec. Ill and IV, and E BC and M B c are 
given by 

jNBjNC f0O dkk 4 u 2 (kR ) 


Ebc(q 2 ) = 
M B c(q 2 ) = 


I2n 2 m 2 Jo (u> B n + x> k )(u. jqn + 
f NB f N c -oo dk k 4 u 2 (kR) 


fl E {q, k ), 
fc), 


(A7) 

(AS) 


127t 2 777 2 7 0 (cubat + U k ) ( UJ C N + Vfc 

where VL B (q,k) and Q,M{q,k) contain the recoil corrections for the intermediate baryons in 
Fig. 1(b) and are given by 


1 r 1 

&E(q, k) = yy— / dx [(E' + m N )(E + m N )] 


1/2 


4mAr 7-i 

1 lq 2 /2 + k-q —q 2 /2-\-k-q 


1 + 


k 2 - q 2 /A 


tt M (q, k) = 


2771tv \ E' + 777tv 

1 C 


E + 777TV 


4777tv 7-1 


1 + 


chr [(A' + m N ){E + m tv)] 1 ^ 2 s m/v 


(E r + m N )(E + ?77tv) 

U(9 2 )j, 

1 1 


Ai(g 2 


(A9) 


+ 


E' + m N E + m N 


Fi(q 2 


k 2 - q 2 /A 


F 2 (q 2 


(A10) 


(E 1 + m N )(E + '/77tv)_ 

Without the recoil of the intermediate baryons (i.e., with k set to zero) the usual CBM 
results are thus recovered, 


il E (q,k) = F,(q 2 ) ~ r)F 2 (q 2 ) = Gf(q 2 ), 
a M (q,k) = Ufa 2 ) + Ufa 2 ) = G { S(q 2 ). 


(All) 
(A 12) 


Now let us discuss the issue of gauge invariance. In the CBM the baryons are assumed 
to be on-mass-shell, thus it only makes sense to discuss current conservation as a weak 
condition for electromagnetic gauge invariance. In a static treatment, current conservation 
holds trivally. With center-of-mass corrections, the electromagnetic form factors are most 
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conveniently calculated in the Breit frame. Since q° = 0 in this frame, current conservation 
is ensured provided that 


9- W|)|J(0)|JV(-|)> = 0. 


(A13) 


For the quark core [Fig. 1(a)], explicit evaluations in both PY and PT projection methods 
guarantee that the matrix element of the spatial component of the current is proportional 
to cr x q, and thus satisfies Eq. (|A13|) . For Fig. 1(b), as shown in the previous paragraph, all 
terms which are odd in k simply vanish after the angular integration over the loop momen¬ 
tum, and the only surviving term is proportional to crxq, therefore this diagram is separately 
gauge invariant. The proof of gauge invariance for Fig. 1(c) is slightly different since it is 
evaluated in the heavy baryon approximation. By expanding the pion held in a plane wave 
and connecting the pion creation/annihilation operators with the CBM hamiltonian and the 
physical baryon states, it is easy to show that [[1^ 


j w (0) oc / dkk k ■ cr x q = k 2 cr x q. 


(A14) 


Thus this current is also transverse with respect to q. Since the total 
is just the sum of the three contributions [Fig. 1(a), 1(b) and 1(c)] 
current conservation, Eq. (|A13|) , is satisfied in this calculation. 


electromagnetic current 
in the CBM, and hence 
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TABLES 

TABLE I. Magnetic moments of the nucleon. The static case refers to the original CBM 


results without center of mass correction, and PY and PT are for two calculations with momentum 
projected wave functions. The experimental values are 2.79 hb and —1.91 hb , respectively. 



proton 

neutron 

R(frn) 

static 

PY 

PT 

static 

PY 

PT 

0.8 

2.49 

2.25 

2.36 

-2.06 

-1.89 

-1.97 

0.9 

2.44 

2.18 

2.30 

-1.96 

-1.78 

-1.86 

1.0 

2.46 

2.18 

2.31 

-1.92 

-1.73 

-1.81 

1.1 

2.53 

2.23 

2.36 

-1.93 

-1.71 

-1.81 
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FIGURES 


FIG. 1. 
pion loop). 


Fig. 1 



B 

(a) (b) (c) 


Diagrams illustrating the various contributions included in this calculation (up to 
The intermediate baryons B and C are restricted to the N and A. 


one 
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FIG. 2. The effect of the center of mass correction and Lorentz contraction for the charge form 
factor of the bare proton. The bag radius is taken to be 1.0 fm. The “static” curve refers to the naive 
MIT cavity approximation, PY and PT stand for Peierls-Yoccoz and Peierls-Thouless projection, 
respectively, and PYL and PTL for the corresponding versions with the Lorentz contraction. 


Fig.2 
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Fig.3 



Q 2 (GeV 2 ) 


FIG. 3. The individual contributions to the proton charge form factor with the bag radius R = 
1.0 fm. The quark part is calculated using the Lorentz contracted PT wave functions. Experimental 
data are taken from Ref. |25||. 
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Fig.4 



Q 2 (GeV 2 ) 


FIG. 4. The individual contributions to the neutron charge form factor. The key is as in Fig. 3, 
except that the experimental data is from Ref. |2(|. As the contribution from Fig. 1(c) is negative, 
we show its magnitude for convenience. 
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Fig.5 



FIG. 5. The proton charge form factor for three different bag radii. Lorentz contracted PT 
wave functions (with ij p = 1) are used in the calculations. Data are the same as in Fig. 3. 
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Fig.6 



Q 2 (GeV 2 ) 


FIG. 6. The neutron charge form factor using Lorentz contracted PT wave functions. Data are 
the same as in Fig. 4. 
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Fig.7 



Q 2 (GeV 2 ) 


FIG. 7. The proton magnetic form factor using Lorentz contracted PT wave functions. Exper¬ 


imental data are from Ref. [25 
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Fig.8 



Q 2 (GeV 2 ) 


FIG. 8. The neutron magnetic form factor using Lorentz contracted PT wave functions. Ex¬ 


perimental data are from Ref. [27,28]. 







